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Summary 

The  concept  of  the  local  heat  transfer  simulation  (LHTS) 
of  the  high  enthalpy  flow  action  on  a vehicle  stagnation 
point  formulated  in  refs.  1,  2 is  based  on  the  requirement 
to  locally  provide  in  a ground  test  the  same  boundary 
layer  on  the  model  at  the  stagnation  point  as  at  the  re- 
entry conditions.  The  present  methodology  of  the 
extrapolation  from  ground  to  flight  consists  of  the  three 
main  parts:  1)  the  recalculation  of  the  test  conditions  to 
flight  parameters,  2)  the  prediction  of  the  test  conditions 
for  the  given  hypersonic  flight  parameters  and  3)  the 
validation  of  the  extrapolation  procedure. 

The  application  and  validation  of  this  concept  are  carried 
out  through  the  comparative  analysis  of  the  computed 
heat  transfer  rates  and  profiles  of  the  gas  temperature  and 
atoms  fractions  within  boundary  layers  near  the  model  in 
the  subsonic  high  enthalpy  air  flow  and  near  the  blunt 
body  in  the  hypersonic  air  flow,  whose  parameters  are 
extrapolated  from  the  inductively  heated  air  plasma 
experiment.  It  is  established  that  the  LHTS  concept 
provides  an  excellent  accuracy  if  air  temperatures  at  the 
outer  edges  of  both  boundary  layers  near  a model  and  a 
body  are  close  to  equilibrium  values,  or  if  the  surface  has 
a high  or  moderate  catalycity. 

The  algorithm  of  determination  of  the  trajectory  point  for 
which  the  complete  local  duplication  of  the  heat  transfer 
is  possible  is  presented.  Another  examples  of  the  LHTS 
concept  applications  are  shown  through  an  analysis  of 
the  capabilities  of  the  IPG-4  plasmatron  for 
thermochemical  simulation  at  the  Mars  probe  entry 
conditions,  and  through  the  prediction  of  the  plasmatron 
subsonic  test  parameters  for  the  duplication  of  the 
stagnation  point  heat  transfer  to  the  Mars  Pathfinder 
aeroshell  at  the  peak-heating  conditions  in  the  Martian 
atmosphere.  The  computed  stagnation  point  heat  flux 
range  for  the  model  in  the  subsonic  dissociated  carbon 
dioxide  flow  in  the  whole  range  of  surface  catalycity  is 
found  in  sufficient  agreement  with  the  viscous-shock 
layer  analysis  carried  out  in  ref.  3 for  the  Mars  Pathfinder 
aeroshell  without  the  surface  ablation. 

1.  Introduction 

In  fact,  only  partial  heat  transfer  simulation  for 

atmospheric  entry  conditions  could  be  achieved  by  using 
high  enthalpy  wind  tunnels  [ref.  4],  The  stagnation  point 
configuration  is  most  important  for  the  TPM  samples 


testing  and  surface  catalycity  rebuilding  [refs.  5-14].  In 
particular,  the  catalytic  properties  of  the  TPM  for  the 
Buran  space  vehicle  were  studied  by  using  the  inductive 
plasmatron  IPG-2  on  cylindrical  models  of  the  30  mm  in 
diameter  with  a flat  face  [ref.  7],  In  subsonic  high- 
enthalpy  flows  of  dissociated  nitrogen,  oxygen  and  air, 
the  effective  probabilities  ymo  of  the  catalytic 
recombination  of  the  N and  O atoms  on  the  ceramic  tile 
coating  and  the  antioxidation  coating  of  the  carbon- 
carbon  material  were  determined  at  the  enthalpy  20-22 
MJ/kg  in  the  pressure  range  0. 1-0.3  atm  and  the  surface 
temperature  up  to  1750  K:  it  was  found  that  for  the  tile 
surface  ymo<3- 10'1 * 3,  for  the  antioxidation  coating  of  the 
C-C  material  ymo  ~7-10"3  [refs.  2,  7],  Those  ground  data 
were  completely  confirmed  by  the  Bor  and  Buran  space 
vehicles  flight  experiments  [refs.  15,  16], 

The  first  data  on  catalytic  efficiencies  of  quartz-based 
materials  and  C-C  materials  have  been  recently  predicted 
from  heat  transfer  tests  in  dissociated  carbon  dioxide 
flows  performed  by  using  the  100-kW  inductive 
plasmatron  IPG-4  in  subsonic  regimes  [refs.  12-14],  The 
question  appears  how  we  can  extrapolate  those  data  for 
entry  conditions  in  the  Martain  atmosphere? 

It  is  almost  obvious  now  that  the  requirements  to 
duplicate  in  tests  only  full-scale  values  of  stagnation 
pressure  and  surface  temperature  [ref.  9]  or  total  enthalpy 
[ref.  10]  are  insufficient  for  the  complete  simulation  of 
the  heat  transfer  and  oxidation  and  we  still  need  to 
answer  the  key  question:  for  what  hypersonic  flight 
conditions  can  we  substantially  use  the  ground  tests  data 
on  the  TPM  catalytic  properties  or  thermo-chemical 
resistance  for  the  prediction  of  the  full-scale  heat  transfer 
or  vehicle  surface  aging? 

For  the  case  of  the  stagnation  point  heat  transfer  (most 
important  from  the  practical  point  of  view)  this  question 
can  be  solved  on  the  basis  of  the  LHTS  concept 
formulated  in  refs.  1,  2.  This  concept  includes  the 
requirements  to  provide  in  high-enthalpy  tests  the  same 
values  of  the  total  enthalpy,  the  stagnation  pressure  and 
the  velocity  gradient  at  the  stagnation  point  of  the  model 
as  in  a hypersonic  flow  around  the  vehicle  at  the  given 
reentry  trajectory  point.  By  means  of  this  rather  simple 
theory  the  ground  test  data  could  be  extrapolated  to  flight 
conditions  if  we  know  flow  fields  around  a model  and  a 
body  (practically  - effective  radii  of  a model  and  a body 
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nose).  Undoubtedly,  an  accuracy  of  such  extrapolation 
should  be  estimated. 

In  this  paper  some  results  of  application  and  validation  of 
this  concept  are  presented.  The  validation  is  carried  out 
through  the  direct  comparison  of  the  numerical  solutions 
of  the  ID  boundary  layer  problem  for  the  high-enthalpy 
experiment  and  the  corresponding  ID  thin  shock  layer 
problem  for  the  extrapolated  flight  parameters.  The  well 
documented  experimental  and  numerical  results  for 
subsonic  high  enthalpy  air  flow  parameters  in  the 
inductive  plasmatron  IPG-2  [ref.  2]  are  used. 

By  means  of  the  LHTS  theory  the  test  parameters  are 
recalculated  to  hypersonic  flight  conditions  and  then  the 
nonequilibrium  boundary  layer  near  the  model  and  the 
shock  layer  near  the  body  along  the  stagnation  line  are 
computed  for  a 5-species  dissociated  air.  Then  stagnation 
point  heat  fluxes,  temperature  and  the  profiles  of  N and 
O atoms  mass  fractions  within  two  boundary  layers  are 
compared. 

It  is  shown  that  the  surface  catalycity  and  the 
displacement  of  air  temperature  from  the  equilibrium 
values  at  the  outer  edges  of  boundary  layers  are  the 
factors  in  the  actual  accuracy  of  the  presented  test-to- 
flight  extrapolation. 

The  algorithm  of  the  determination  of  the  trajectory  point 
for  which  the  local  thermochemical  simulation  could  be 
realized  by  using  plasmatron  is  described.  Two  another 
examples  of  the  applications  of  the  LHTS  concept  are 
demonstrated  through  analysis  of  the  Mars  probe 
trajectory  [ref.  17]  and  test  parameters  required  for  the 
stagnation  point  heat  transfer  simulation  for  the  Mars 
probe  and  the  Mars  Pathfinder  aeroshell  at  the  trajectory 
peak-heating  point  in  the  Martian  atmosphere  [ref.  3]  by 
using  the  100-kW  inductive  plasmatron  IPG-4  [refs.  11- 
14]  and  the  standard  European  model  configuration  [refs. 
9,  12-14].  The  numerically  predicted  heat  transfer  range 
for  the  experiment  in  a subsonic  high-enthalpy  carbon 
dioxide  flow  is  found  in  good  agreement  with 
computations  for  the  Mars  Pathfinder  aeroshell  from  ref. 
3 in  the  whole  range  of  the  surface  catalycity. 

2.  Concept  of  the  Local  Heat  Transfer  Simulation 

Our  way  to  study  the  problem  of  the  ground-to-flight 
extrapolation  lies  in  the  use  of  the  theory  of  the  local  heat 
transfer  simulation  formulated  in  refs.  1,  2.  At  least  for 
the  case  of  the  stagnation  point  heat  transfer  we  can  point 
out  the  hypersonic  flow  parameters  and  a blunt  body 
radius  which  are  in  direct  correspondence  with  high 
enthalpy  subsonic  flow  parameters  and  a model  radius 
(Fig.  1 ) if  both  surfaces  have  the  same  catalycity  and 
emissivity  and  also  the  same  heat  transfer  boundary 
conditions  (for  example,  radiative-equilibrium  walls). 

The  analysis  is  based  on  the  boundary  layer  theory  for 
dissociated  reacting  gases.  Fay  and  Riddell  theory  [ref. 
1 8]  gives  the  following  expression  for  the  heat  flux  at  the 


Fig.  1.  Schemes  of  hypersonic  flow  around  a blunt  body  and 
the  stagnation  point  heat  transfer  simulation  in  subsonic  high 
enthalpy  flow  in  plasmatron. 

stagnation  point  of  a blunt  body  with  a fully  catalytic 
surface  which  is  being  flown  by  a hypersonic  flow: 


% = 0.763  Pr~0-6 

*(ht  - V 


ds 


pa 
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(2-1) 


For  the  case  of  the  frozen  boundary  layer  the  Goulard’s 
solution  of  the  boundary  layer  problem  [ref.  19]  gives  the 
following  formula  for  the  stagnation  point  heat  flux  to 
the  wall  with  arbitrary  catalycity: 


1w  = ‘(H,  - V[l  + 

+(Le2/3<p  - 1 ]. 


\\/2,dJJe  ,1/2/ 
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1 + 


A-1 


Here  H is  the  total  enthalpy,  h is  the  enthalpy,  h°A  is  the 

enthalpy  of  formation  of  atoms,  c is  the  mass  fraction  of 
atoms,  p is  the  density,  p is  the  viscosity,  (dU/ds)0  is  the 
velocity  gradient,  Pr  is  the  Prandtl  number,  Le  the  Lewis 
number,  Sm  is  the  Schmidt  number.  The  subscripts  e and 
w denote  the  outer  edge  of  the  boundary  layer  and  the 
body  surface,  respectively. 


The  analysis  of  the  formulae  (2.1)  and  (2.2)  has  revealed 
the  three  independent  factors  which  control  the  heat  flux: 
He-h  w,  pt  and  (dUJds)0.  In  the  case  when  the  flow  at  the 
outer  edge  of  the  boundary  layer  is  in  equilibrium,  pe  is 
the  function  of  he(=He)  and  pe,  and  accordingly  to 
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Goulard’s  solution  [ref.  19]  /z^ appears  to  be  a parameter 
dependent  on  he,  pe,  (dUJds)0  and  Tw.  Therefore,  the 
stagnation  point  heat  flux  qw  is  completely  determined  by 
the  boundary  condition  on  the  body  and  the  three 
parameters  at  the  outer  edge  of  the  boundary  layer:  the 
enthalpy,  the  stagnation  pressure  and  the  velocity 
gradient. 

We  see,  that  the  sufficient  conditions  of  equality  of  the 
heat  fluxes  in  flight  and  in  ground  test  are  the  equalities 
of  these  three  parameters  in  the  two  flows: 

Ms  Hoc,  Pcs  P v x Piv-.r  (dUJds) os  (dUe/ds)ox  (2.3) 

where  the  subscripts  x and  S relate  to  the  flight  and 
ground  conditions. 

Let’s  consider  the  conditions  of  duplication  of  the 
stagnation  point  heat  transfer  to  a blunted  body  of  radius 
Rw,  which  is  being  flown  by  a hypersonic  flow  with 
velocity  V„  and  density  pm<  in  an  axisymmetric  high 
enthalpy  subsonic  flow  with  velocity  Vs  and  static 
pressure  ps  by  using  a cylindrical  blunted  model  of 
radius  Rn,  {Fig.  1). 

From  the  first  equality  (2.3)  the  equality  of  the  total 
enthalpies  of  the  two  flows  follows: 

<2-4> 


The  second  necessary  condition  reduces  to  the  relation 
between  the  static  pressure  of  the  free  stream  in  a plasma 
wind  tunnel  and  the  given  stagnation  pressure  pw,  which 
is  described  by  the  approximate  Poisson’s  adiabat 
equation  for  a real  gas  [ref.  20]: 


ps(\  + M2s/*S/(7*S~l)  = pw , (2.5) 

1 

J*s  1-VpA 


where  y,s  is  the  effective  specific  heat  ratio,  and  Ms<\  is 
the  Mach  number  of  the  simulating  flow. 

Let’s  represent  the  third  condition  - the  equality  of  the 
stagnation  point  gradients  - in  the  form 


V„  / R 

S m 


V / R*, 

oc  N’ 


(2.6) 


R = V,.  / (dU  / ds)n„ 
m,N  S.cc  ' e /0<S,cc 


(2.4)-(2.5)  we  obtain  the  parameters  of  the  ground  free 
stream  in  the  following  form: 


vs/vx^,  hs/Hx= i-e, 

ps/pw=a-?/*s/(y*s-v, 


M>=—? 

^R*m/R*N,^(vl/2Hje, 


'*S 


' *s 


+ 1 


(2.7) 

(2.8) 

(2.9) 

(2.10) 


The  expressions  for  calculating  the  velocity,  enthalpy, 
pressure  and  Mach  number  of  the  simulating  flow  (2.7) 
are  universal  in  form:  they  contain  only  two 

dimensionless  parameters  £ and  C,.  For  the  hypersonic 
flight  in  an  atmosphere  F2  / 2Hx  = 1 , then  £=£. 


Thus,  for  high-enthalpy  flows  in  which  the  stagnation 
point  heat  fluxes  at  the  body  and  the  model  are  equal,  the 
ratios  V„  / Vx  and  h/Hx  are  completely  determined  by 

the  elementary  universal  relations  (2.7).  The  relations 
(2.8)  and  (2.9)  for  high-enthalpy  flows  are  the  universal 

implicit  relations  of  the  functions  pi/Q  = PS(Q/  Pw 

and  MS(Q  , since  y,s  depends  on  the  values  of  ps  and 

hs..  In  Fig.  2 we  have  plotted  the  functions  pidQ  and 

MS(Q  calculated  for  equilibrium  air  at  7/^=3 2 (a),  16 

{ b ) and  8 (c)  MJ/kg.  The  curves  1-3  correspond  to  pw 
=10"3,  1 and  102  atm.  The  broken  curves  correspond  to  a 
perfect  gas  with  constant  y=l.4. 


where  R*  and  R*N  are  the  effective  radii  of  the  model 
and  the  body  at  the  stagnation  points.  Assuming  the 
parameters  Vx,  Hx,  pw,  R*m , and  R*N  to  be  known  from 


Fig.  2.  Universal  dependencies  of  the  dimensionless  static 
pressure  and  Mach  number  of  simulating  subsonic  air  flow 
versus  similarity  parameter  f [ref.  1], 

As  may  be  seen  from  Fig.  2,  the  real  properties  of  highly 
dissociated  air  essentially  affect  the  functions 
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p$(Q  and  M^CQ , which  are  quite  important  for  the 
analysis  of  the  test  conditions. 
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f \ 1/2 
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3.  Recalculation  of  Test  Parameters  to  Flight 
Conditions 

For  the  experimental  practice  of  the  heat  transfer 
simulation  the  case  £«\  is  very  important:  practically, 
this  case  occurs  when  Rm«RN.  Let’s  consider  this  case  in 
detail  taking  into  consideration  that  for  a hypersonic 
axially  symmetric  flow  around  a smooth  blunted  body  we 
have  the  following  well-known  formulae  (ref.  20): 


H = —V2,  p„/  = P V2.  (3-1) 

oc  ^ 00  r°°  00  V 7 


1 


r*,i 

1 Vn_ 

l ds  j 

0 

v3  P e) 

Rn 

Then  for  the  ground  subsonic  free  stream  conditions  we 
have: 


Ps-vJl 

(3.3) 

fs  p y/2  r‘ 

u r cc  m y 

^3  p J Rn 

(3.4) 

In  the  case  of  subsonic  jet  the  subscript  S denotes  the 
center  of  the  flow. 


As  an  example,  let’s  analyze  the  well  characterized 
subsonic  regime  performed  in  the  IPG-2  plasmatron  for 
the  air  flow  at  the  pressure  0.1  atm  , the  generator  anode 
power  37.4  kW  and  the  mass  flow  rate  through  the 
discharge  channel  2.8g/s  [ref.  2],  At  the  distance  of 
30  mm  from  the  plasmatron  exit  section,  the  enthalpy 
T/e=2.19T07  m2/s2 , the  velocity  F/M80  m/s,  the  density 
p=  3.86  10'3  kg/m3,  Reynolds  number  Res=psVsRJps= 
s58.6,  Mach  number  A/s-=  0.14.  The  radius  of  the  IPG-2 
plasmatron  channel  Rc=3  10'2  m and  for  the  cylindrical 
model  of  the  radius  Rm=  1.5-1  O'2  m we  have 

/?^=1.27?m=1.8T0"2  m,  in  accordance  with  (3.5). 

For  these  subsonic  air  flow  parameters  we  have 
determined  the  parameters  of  the  hypersonic  air  flow  and 
the  nose  radius  of  a blunt  body  from  expressions  (3.6): 
F„=6 620  m/s,  /v=2.28T0'4  kg/m3,  Rp=  0.265  m.  The 
density  value  corresponds  to  the  altitude  62.4  km  in  the 
Earth  atmosphere.  Correspondingly,  the  flight  Reynolds 
number  ReK=l75 0 and  the  Mach  number  Ma=20.  We 
see,  that  Reynolds  and  Mach  numbers  are  quite  different 
for  those  ground  test  and  hypothetical  flight  regime.  That 
means  that  Reynolds  and  Mach  numbers  are  not  the 
similarity  parameters  for  the  stagnation  point  heat 
transfer. 


The  conditions  (3.3)  are  rather  simple,  but  the  condition 
(3.4)  is  not  trivial  and  contains  the  geometrical  parameter 

- effective  radius  R*m,  which  depends  on  the  channel 

radius  Rc,  the  model  radius  Rm  and  should  be  computed 
for  the  test  configuration  by  using  CFD  methods.  For 
subsonic  jets  over  cylindrical  models  with  a flat  face  we 

obtained  the  following  approximation  for  R*m  from 

numerical  solutions  of  the  Navier-Stokes  equations 
[ref.  1] 


r: 


* 


m 


m 


2-1- 1,68(1  - 1)2  - 1,28(1  -l)3 

Rc.  l>l  0 = Rm/Rc)  ~ 


1<1 


(3.5) 


In  accordance  with  (3.3),  (3.4)  there  is  one-to-one 
correspondence  between  subsonic  jet  parameters 
hs,ps(~  pw),Vs  and  the  given  dimensions  Rc  and  R„„ 

on  the  one  hand,  and  the  parameters  of  the  hypersonic 
flow  and  the  nose  radius  of  a blunt  body,  on  the  other 


V 

OO 


(3.6) 


4.  Method  of  the  LHTS  Concept  Validation 

The  proposed  method  of  the  LHTS  concept  validation 
consists  in  the  direct  comparison  of  the  calculated  heat 
transfer  rates  and  the  profiles  of  the  temperature  and 
species  fractions  within  the  boundary  layers  near  the 
stagnation  points  of  the  model  and  the  body  for  subsonic 
and  hypersonic  flows  conditions  linked  by  the  correlation 
(3.6).  Such  a comparison  is  carried  out  below. 

A method  for  the  calculation  of  the  heat  transfer  rates  at 
the  stagnation  point  of  a model  with  a flat  face  exposed 
to  a subsonic  jet  of  a viscous  multicomponent  reacting 
gas  was  developed  in  refs.  2,  7,  21,  22,  13.  Here  we  also 
use  the  concept  of  a boundary  layer  with  finite  thickness 
in  the  vicinity  of  the  stagnation  point.  The  thickness  8 of 
the  boundary  layer  and  the  hydrodynamic  parameters  at 
its  outer  edge,  including  the  flow  vorticity,  which  are 
given  below  in  the  formulation  of  the  problem,  are 
determined  from  the  numerical  solution  of  the  full 
Navier-Stokes  equations,  which  is  considered  as  an  outer 
solution. 

In  the  case  of  a one-temperature  multicomponent  mixture 
of  atoms  and  molecules  in  the  equilibrium  excitation  of 
the  vibrational  degrees  of  freedom,  the  system  of  the 
ordinary  differential  equations  describing  the  flow  within 
the  boundary  layer  near  the  stagnation  point  of  the  model 
has  the  following  form  (ref.  21): 
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A = 8 / 7?  , u = u,  / u,  , u,  = du / dx,  u = -F°  / V , 
m 1 \e  1 s 


a = -v  (du,  / dy)  / uf,  U = U°/V,  x = x°  / R , 
e e 1 7 e le  s m 


For  our  case  A=8/Rm=0A,  oc=0.52,  M|c,  =0.445,  ae= 2.626. 
It  was  assumed  in  calculation  that  Pr=0.71,  Sc= 0.65  and 
ju~T01\ 

We  assume  that  the  following  gas-phase  reactions  in  a 5- 
species  air  occur  within  the  boundary  layer:  1)  02+M+-> 
O+O+M,  2)  N2 + M +N+  M,  3)  N0+M-f->N+0+M,  4) 
0+N2<->N+N0,  5)  0+N0<->N+02.  The  chemical 

equilibrium  is  supposed  to  be  at  the  outer  edge  of  the 
boundary  layer:  />y=0.1  atm,  He= 2.19-107  m2/s2, 

re=5960K,  Cjv2=0.4576,  CO2=0.3513-10’3,  CN2= 

0.2922-1  O'2,  Cn=  0.3091,  Co=0.2304. 

The  thin  viscous  shock  layer  model  was  used 
independently  for  the  computation  of  the  hypersonic  air 
flow  (VrJ~6620  m/s,  /7M=2,28-10'4  kg/m3)  past  a sphere  of 
the  radius  P^0,265  m.  At  the  formulation  of  this 
problem  the  same  boundary  conditions  at  the  wall  and 
the  same  chemical  reactions  rates,  as  for  the  subsonic 
flow,  were  used.  For  the  numerical  solutions  of  the 
problem  (4.1),  (4.2)  and  the  viscous  shock  layer  problem 
the  fourth-order-accurate  finite-difference  scheme  was 
used.  All  computations  were  made  for  the  surface 
temperature  7^=1500  K. 


y = y°/R,  p = p°/p,  T=T°/T,  h.  = h?  / H . 

77  m’  r r re  e’  i i e 

The  boundary  conditions  at  the  outer  edge  of  the 
boundary  layer  and  on  the  surface  of  the  model  are: 

n = 1:  u=H=  1,  c.  = c (i  = - N ) (4.2) 

1 ‘e  e 


5.  Results  of  Validation  for  Subsonic  Air  Test 

Fig.  3 shows  the  calculated  dependencies  of  the 
stagnation  point  heat  fluxes  qw  to  the  model  (M«  1)  and 
the  body  (M»  1 ) as  the  functions  of  the  heterogeneous 
recombination  rate  constant  Kw  for  the  case  KwctK-wn 
=KW  at  test  and  flight  conditions  specified  above. 


tl  = 0:  « = / = 0,  T=Tw,  y = 0, 


K = —22-  — x- 
wi  2-y.pnm.  J 

In  (4.1),  (4.2)  If,  V°  are  the  velocity  components  in  the 
cylindrical  coordinate  system  x°,  y°  which  is  related  with 
a flat  face,  p is  the  density,  c,  is  the  mass  fraction,  h,  is 
the  enthalpy,  h,  is  the  energy  of  formation  of  the  species 
/;  c' , J*  are  the  mass  fraction  and  the  diffusive  flux  of 
the  chemical  element  j;  H is  the  enthalpy  of  the  gas 
mixture,  T is  the  temperature,  m is  the  molecular  weight, 
KWi  and  y,  are  the  effective  heterogeneous  recombination 
rate  constant  and  the  catalytic  efficiency;  k is  the 
Boltzmann  constant,  m,  is  the  molecular  weight  of  the 
species  i,  N is  the  number  of  species,  Ne  is  the  number  of 
chemical  elements,  and  77  is  the  Dorodnitsyn’s  variable. 
The  circle  superscript  denotes  dimensional  quantities. 

In  the  momentum  equation  the  parameter  ae=const  takes 
into  account  the  vorticity  of  the  flow  at  the  outer  edge  of 
a boundary  layer  of  the  thickness  8.  Parameters  ae,  Ue 
and  uu  are  determined  from  the  profiles  of  the  velocity 
components  obtained  from  the  numerical  solution  of  the 
Navier-Stokes  problem  for  a viscous  reacting  gas  jet  flow 
past  a cylinder  with  the  flat  face  at  M«  1 [refs.  2,  7,  19]. 


K 


\Reu, 


ScA 


-C  = 
/ 


VsUU  V' 


Kw 

Fig.  3.  Stagnation  point  heat  transfer  rates  qw,  (W/cm2)  versus 
effective  rate  constant  of  surface  atom  recombination 
Kfy,(m/s)  for  plasmatron  test  (M«l)  and  extrapolated  to  flight 
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For  a highly  catalytic  surface  the  difference  in  the  heat 
fluxes  qw  under  corresponding  conditions  for  the 
subsonic  and  the  supersonic  flows  is  only  about  5%.  As 
Kw  decreases,  this  difference  increases  and  for  a 
noncatalytic  surface  in  the  subsonic  jet  the  heat  flux  is 
30%  less  then  the  heat  flux  in  the  corresponding 
hypersonic  flow. 

So,  we  see,  that  the  accuracy  of  the  heat  transfer 
simulation  on  the  basis  of  the  LHTS  theory  developed  in 
refs.  1,  2 depends  on  the  surface  catalytic  efficiency.  The 
accuracy  is  rather  good  for  surfaces  with  high  and 
moderate  catalycity  but  it  looks  insufficient  for  a 
noncatalytic  one.  Nevertheless,  both  curves  in  Fig.  3 are 
functionally  similar  and  they  have  the  same  practically 
important  interval  10'1  <Kw<  102  m/s,  where  the  heat 
transfer  rates  drastically  depend  on  wall  catalycity. 

Moreover,  it  is  easy  to  find  by  using  Fig.  3 that  at  least 
one  function  exists,  which  is  duplicated  with  a quite 
perfect  accuracy  - the  normalized  heat  flux 


% 


‘V  qWn 
Qwk  _ qWn 


(5.1) 


where  qm  is  the  heat  flux  to  a fully  catalytic  wall  and 


qWn  is  the  heat  flux  to  a noncatalytic  wall. 


(see  Fig.  5):  the  mass  fraction  cN  is  significantly  less  then 
the  equilibrium  value  in  the  subsonic  flow.  On  the  other 
hand,  close  to  fully  catalytic  surfaces  (Fig.  5),  the 
profiles  of  the  N atoms  fractions  for  the  considered 
subsonic  and  hypersonic  air  flows  are  quite  similar 


Y,  cm 


From  here  we  can  conclude  that  the  data  on  the  effective 
catalytic  rates  for  atomic  oxygen  and  nitrogen 
recombination  on  the  Buran  TPM  \<KW<3  m/s,  obtained 
in  subsonic  jets  [refs.  2,  7,  20],  are  quite  applicable  to 
hypersonic  re-entry  conditions  in  the  Earth  atmosphere 
with  flight  parameters  of  the  same  orders  of  magnitude 
that  calculated  above. 

For  the  more  clear  understanding  of  the  LHTS 
capabilities  we  will  compare  the  profiles  of  the  air 
temperature,  N and  O atoms  fractions  across  the 
boundary  layer  near  the  model  and  across  the  shock  layer 
near  the  body  at  the  same  conditions  for  two  limiting 
cases:  fully  and  noncatalytic  walls. 

Fig.  4 shows  the  temperature  profiles  within  the 
boundary  layer  near  the  model  and  within  the  shock  layer 
near  the  blunt  body  for  a fully  catalytic  wall  case.  We  see 
that  temperature  profiles  are  rather  different  at  the  outer 
edges  of  two  boundary  layers.  For  the  hypersonic  flow 
the  temperature  Te  considerably  (~  by  2000  K)  exceeds 
the  equilibrium  value  because  the  shock  layer  is  fully 
nonequilibrium  in  this  case,  but  at  the  same  time  the  two 
temperature  profiles  are  quite  close  to  each  other  near  the 
walls.  This  leads  to  good  agreement  between  thermal 
conductive  parts  of  the  heat  fluxes  to  fully  catalytic 
surfaces  of  the  body  and  the  model. 


Fig.  4.  Temperature  profiles  along  stagnation  line  (fully 
catalytic  wall):  I - the  boundary  layer  near  the  model  in 
subsonic  flow  (Rm-1.5  10~2  m);  2 - the  hypersonic  shock  layer 
near  the  body  (R^O.265  m). 


N 


0.5  Y cm  1.0  1.5 


Within  the  hypersonic  shock  layer  we  can  see  some 
displacement  of  the  N atoms  mass  fraction  from  the 
equilibrium  value  at  the  outer  edge  of  the  boundary  layer 


Fig.  5.  Profiles  of  the  mass  fractions  of  N and  O atoms  along 
stagnation  line  (fully  catalytic  wall):  1 - the  boundary  layer 
near  the  model  in  subsonic  flow,  2 - the  hypersonic  shock  layer 
near  the  body. 


8B-7 


including  the  fine  details  of  different  behavior  of  the 
profiles  of  N and  O atoms  due  to  the  exchange  reactions 
in  dissociated  air  and  mentioned  above  in  the  chapter  4. 
In  fact,  the  contributions  in  the  heat  fluxes  due  to 
nitrogen  atoms  diffusion  are  insignificant  in  this  case. 
Also  we  can  see  on  Fig.  5 that  the  excellent  simulation 
accuracy  is  achieved  for  the  atomic  oxygen  fraction 
profile  within  the  whole  boundary  layer  because  oxygen 
is  completely  dissociated  at  the  outer  edges  of  both 
boundary  layers.  As  a result  the  contributions  in  the  heat 
fluxes  due  to  the  atomic  oxygen  diffusion  for  both  flows 
are  equal  within  5%.  The  same  is  correct  for  the  total 
heat  fluxes  to  a fully  catalytic  wall. 

For  the  noncatalytic  wall  case  the  situation  is  more 
dramatic  as  we  can  see  in  Fig.  6 and  7.  The  temperature 
profiles  across  boundary  layers  are  different  exterior  to 
the  nearest  vicinities  of  stagnation  points  and  the  heat 
flux  caused  by  the  to  thermal  conductivity  is  30  % higher 
in  the  supersonic  flow,  then  in  subsonic  one. 


Fig.  6.  Temperature  profiles  along  stagnation  line  (non- 
catalytic wall):  1 - the  boundary  layer  near  the  model  in 
subsonic  flow,  2 - the  hypersonic  shock  layer  near  the  body. 

Both  boundary  layers  are  almost  frozen  and  diffusion 
does  not  influence  heat  transfer.  The  atomic  nitrogen 
fraction  near  the  surface  of  the  model  is  higher  in  the 
plasmatron  test,  but  the  atomic  oxygen  fraction  profiles 
within  boundary  layers  are  quite  close  in  test  and 
hypothetical  flight  (Fig.  7).  That  means  the  formulae 
(3.6)  for  ground-to-flight  extrapolation  should  provide  a 
complete  simulation  of  the  diffiision  flux  and  the  partial 
pressure  of  atomic  oxygen  and,  therefore,  surface 
catalysis  and  oxidation  processes.  These  are  the 
arguments  for  using  here  the  term  “thermochemical” 
simulation. 


Fig.  7.  Profiles  of  the  mass  fractions  of  N and  O atoms  along 
stagnation  line  (noncatalytic  wall):  1 - the  boundary  layer 
near  the  model  in  subsonic  flow,  2 - the  hypersonic  shock  layer 
near  the  body. 

It  was  established  in  ref.  23,  that  an  accuracy  of  the  heat 
transfer  duplication  for  low  catalytic  surface  was 
improving  when  the  pressure  in  subsonic  high  enthalpy 
flow  was  increasing.  When  ps>0.2  atm  the  difference  in 
the  qw  for  subsonic  and  hypersonic  flow  was  not  more 
5%  for  surfaces  with  AV>1  m/s  (the  quite  practical  case). 
The  example  of  the  excellent  duplication  of  the 
temperature  distribution  across  the  boundary  layer  for 
hypersonic  flow  conditions  we  can  see  on  Fig.  8 for  high 
enthalpy  subsonic  test  in  dissociated  nitrogen  [ref.  23], 


Fig.  8.  Temperature  distributions  along  stagnation  line  in 
hypersonic  shock  layer  (1:  Mx=20,  Vx=6490  m/s,  Z=62  km, 
Rw=0.95  m)  and  in  subsonic  jet  (2:  Ms=0.04,  Vs=50m/s, 
Ps=0.3  atm,  Fls=21  MJ/kg,  Rm=1.5-10~2m). 
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We  can  conclude  that  the  displacement  of  air 
temperature  from  the  equilibrium  value  at  the  outer  edges 
of  boundary  layers  in  high-enthalpy  tests  or  in  flight  is  a 
factor  in  actual  accuracy  of  the  LHTS  concept.  One 
should  expect  that  the  agreement  between  the  heat  fluxes 
and  two  boundary  layer  structures  will  be  improved  when 
the  air  temperatures  outside  boundary  layers  are  closer  to 
equilibrium,  for  example,  - for  higher  values  of  the 
pressure  as  has  already  been  predicted  for  the  high- 
enthalpy  nitrogen  experiment  in  ref.  23 . 

6.  The  Trajectory  Point  for  the  Complete  Local 
Duplication  of  Heat  Transfer 

One  in  the  main  features  of  the  subsonic  high  enthalpy 
jets  is  the  nonuniformity  of  the  enthalpy  and  velocity 
profiles  at  the  plasma  generator  channel  exit  and  the 
decreasing  of  these  characteristics  along  the  flow  axis.  In 
order  to  take  these  effects  into  account,  in  relations  (2.7) 
instead  of  Vs  and  hs  we  should,  substitute  the  velocity  and 
enthalpy  values  on  the  axis  of  the  free  stream 
corresponding  to  the  distance  L from  the  plasmatron 
channel  exit  to  the  model. 

For  the  modeling  in  a hyposonic  flow  (M«  1),  taking 
into  account  the  above  we  can  generalize  relations  (3.3), 
(3.4)  in  the  following  form 

hs  - i, cJ(L)rl  ps  . a-K)p  Jl K - (6.1) 

y e 

, (s  y/2  r* 

<62) 

Here,  the  factors  <pH(L)  and  <s?y( L)  take  into  account 

the  decreasing  of  the  enthalpy  and  velocity  along  the  axis 
of  the  free  subsonic  stream,  the  subscript  S relates  to  the 
center  of  the  plasmatron  channel  exit. 

As  we  mentioned  above,  the  parameters  hs,  ps  and  Vs  are 
in  functional  coupling.  For  a wide  range  of  the  subsonic 
tests  conditions  in  the  optimum  discharge  burning  regime 
these  parameters  could  be  presented  with  a functional 
relationship 


Psvs  = *(hs’&ps)  (6-3) 


where  Q is  the  mass  flow  rate  and  %(hs,Q,ps)  is  the 

specific  functional  characteristic  of  the  plasmatron.  This 
characteristic  should  be  determined  for  each  facility 
which  is  used  for  the  heat  transfer  duplication. 

If  function  x is  known,  we  can  eliminate  ps  and  Vs  in 
(6.3)  using  (6.1)  and  (6.2)  and  obtain  the  following 
relationship  between  the  hypersonic  flow  parameters  p,x 
and  Vx 


p«  - (~‘r',2a -KrVpi)^3  x 


XXC — ,&(T-K)p  V1) 

2qH(L) 


(6.4) 


For  the  following  step  of  the  procedure,  it  is  necessary 
to  present  the  trajectory  of  the  body  in  the  atmosphere  in 
the  parametric  form  Vx-F(px).  The  intersection  of  this 
curve  with  the  curve  (6.4)  in  the  plane  px-Vx,  if  it  exists, 
gives  us  the  trajectory  point  for  which  the  complete  local 
heat  transfer  duplication  could  be  achieved. 

For  the  specified  parameters  Vx,  px,  and  RN.  the 
complete  local  simulation  of  the  heat  transfer  is  possible 
for  the  appropriate  model  (or  channel)  dimension.  In  this 
case  the  test  conditions  can  be  predicted  in  accordance 
with  the  following  algorithm.  The  free  stream  parameters 
hs  , ps,  and  Vs  are  calculated  from  (6.1)  and  (6.3),  and 

then  the  effective  radius  R*  is  calculated  from  (6.2).  The 

necessary  value  of  R*m  is  insured  by  choosing  a model  or 

channel  with  the  appropriate  geometry  (for  models  with 
flat  face  in  accordance  with  expression  (3.5)). 

7.  Analysis  of  the  Mars  probe  trajectory  and 
Requirements  for  Plasmatron  Tests 

Let’s  apply  the  LHTS  concept  in  order  to  estimate  the 
IPG-4  plasmatron  capabilities  for  complete  local 
duplication  of  the  stagnation  point  heat  transfer  for  the 
Mars  probe  trajectory  parameters  in  a subsonic  test  with 
carbon  dioxide.  The  IPG-4  operating  envelope  for  the 
subsonic  regime  with  carbon  dioxide  as  a working  gas 
(Fig.  9)  contains  the  peak-heating  parts  of  the  Mars 
probe  trajectory  [ref.  17], 


Fig.  9.  The  IPG-4  operating  envelope  in  the  stagnation 
pressure-enthalpy  coordinates  for  subsonic  regime  with 
carbon  dioxide,  and  the  Mars  Pathfinder  and  the  Mars  Probe 
trajectories;  the  stars  indicate  the  peak-heating  points  [ref.  24], 
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The  velocity  Vm  the  altitude  H and  the  calculated 
convective  heat  flux  to  fully  catalytic  radiative- 
equilibrium  wall  qw  at  the  stagnation  point  for  the  probe 
with  T^O.8  m are  presented  in  Fig.  10  as  functions  of 
the  entry  time  [ref.  17]. 


Fig.  10.  Mars  probe  trajectory  and  convective  heat  flux  at 
stagnation  point  to  fully  catalytic  radiative-equilibrium  wall 
(R  11=0.8  m.  £,h=0. 85)  [ref.  24]:  1 - H-100,  2 -Vm,  3 - qw. 

In  accordance  with  these  data  and  the  density  distribution 
in  the  Martian  atmosphere  from  ref.  17  we  have 
calculated,  by  using  (3.3),  the  dependencies  p/t)  and 
h/t)  that  are  presented  in  Fig.  77  by  the  curves  7 and  2. 
Also,  in  Fig.  77  the  part  of  the  IPG-4  operating  envelope 
in  coordinates  ps  -hs  for  the  carbon  dioxide  gas  is  shown. 


Fig.  11.  Heat  transfer  parameters  for  the  Mars  probe  within 
the  IPG-4  operating  envelope  as  the  functions  of  the  entry  time 
t (sec)  [ref.  24],  Curves:  1 - enthalpy  (MJ/kg),  2 - normalized 
stagnation  pressure  (pO=l  atm),  3,  4 - velocity  gradients 
(I03  s'1)  in  hypersonic  flow  for  RN=0.8  and  1.25  m,  5 - velocity 
gradient  (103  s'1)  in  subsonic  jet  for  the  euromodel;  points:  a - 
the  peak  heating  point,  b and  c -the  points  of  the  complete 
thermochemical  simulation. 


We  see,  that  the  IPG-4  envelope  contains  the  both  curves 
ps(t)  and  hs(t)  within  the  entry  time  interval  25</<44  s. 
This  interval  includes  the  heat-intensive  part  of  the  probe 
trajectory  during  the  period  30</<44  s (Fig.  10).  The 
maximum  value  of  a heat  flux  should  be  achieved  at  the 
point  t*=38.5  s where  Z=43.0  km  and  ^*=5995  m/s.  One 
could  expect  that  for  an  accurate  simulation  the  heat 
transfer  at  that  most  important  trajectory  point  it  is 
necessary  to  have  the  subsonic  flow  in  the  IPG-4  with  the 
enthalpy  hs  =14.02  MJ/kg  and  the  static  pressure 
/>.v*=6.7-10~2atm. 

Taking  into  account,  that  in  the  IPG-4  subsonic  regimes 
the  pressure  ps  and  the  enthalpy  hs  could  be  controlled 
independently  by  using  the  vacuum  pump  system  and  by 
the  variation  of  energy  input  in  plasma  [refs.  12-14], 
from  Fig.  77  we  can  conclude  that  it  is  possible  to 
duplicate  simultaneously  both  functions  ps(t)  and  h/t)  in 
the  real  time  scale  during  19  seconds. 

But  in  fact  this  way  of  the  duplication  of  only  two  heat 
transfer  parameters  would  not  be  the  complete  simulation 
of  heat  transfer  because  the  satisfaction  of  conditions 

(3.3)  does  not  guarantee  the  satisfaction  of  the  condition 

(3.4) .  Let  us  emphasize  that  two  parameters  of  the  heat 
transfer  simulation  presented  on  Fig.  77  by  the  curves  7 
and  2 depend  neither  on  a body  shape  nor  on  a model 
geometry.  On  the  contrary,  the  equality  of  the  two 
velocity  gradients  VfRm*=VJRN*  depends  on  subsonic 
and  hypersonic  flows  geometry,  shapes  and  dimensions 
of  a body  and  a model. 

In  accordance  with  Fig.  10,  along  the  Mars  probe 
trajectory  within  the  interval  25</<44  s the  velocity 
gradient  at  the  shield  stagnation  point  $ej=Vaft)/RN* 
decreases  monotonously  in  the  range  4.4-103-3.2-103  s'1 
for  Rf=0.S  m and  in  the  range  2.8-103-2.1-103  s"1  for 
77^=1.25  m (see  the  curves  3 and  4 on  Fig.  77). 

Let’s  consider  now  the  variation  of  the  velocity  gradient 
P ep~Vs(t)/Rm*  at  the  stagnation  point  of  the  euromodel 
(the  cylindrical  model  of  50  mm  in  diameter  with  a flat 
face  and  the  rounded  edge  of  11  mm  in  radius  with 
Rm  =3.2-10'2m)  if  the  parameters  ps  and  hs  change  in  the 
way,  as  Fig.  11  requires. 

It  is  very  important  to  understand  that  in  the  plasmatron 
flow  velocity  Vs  is  not  an  independent  parameter  - it  is 
definitely  linked  with  the  gas  flow  rate  through  the 
discharge  channel  Q,  the  enthalpy  hs  and  the  pressure  ps. 
For  further  estimations  we  need  to  have  some 
approximation  for  the  dependence  Vs=Vs(Q,hs,pJ.  As  a 
next  step  on  the  basis  of  the  previous  experimental  and 
numerical  data  we  assume  the  following  approximation 
at  pressures  ps<0. 1 atm 


s 
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The  experimental  dependency  Vs°  (hs)  at  the  pressure 
0.1  atm  is  taken  from  refs.  11,  12. 

For  the  cylindrical  model  of  50  mm  in  diameter  with  the 
rounded  edge  of  1 1 mm  in  radius  (the  euromodel)  and 
the  subsonic  flow  configuration  in  the  IPG-4  (7?c=40  mm) 
on  the  basis  of  the  previous  numerical  solution  of  the 
Navier-Stokes  equations  [ref.  12]  we  have  the  following 
approximations  for  the  velocity  gradient  and  effective 
radius 


dU_ 


ep 


dr 


' 0 


V 

= 0.78—^, 
^ m 


R*  = 1.287?  =3.2-10  2 m 
m m 


(7.2) 


In  Fig.  11  the  curve  5 shows  the  velocity  gradient 
P ep(t)~Vs(t)/Rm*  that  has  been  calculated  from  formula 
(7.2)  taking  into  account  (7.1)  and  dependencies  ps(t) and 
h/t).  We  can  see,  that  this  velocity  gradient  p ep(t) 
decreases  monotonously  within  the  time  interval 
35<t<44  s from  4.7-103  to  1.15T03  s'1.  The  curve  5 
crosses  the  curves  3 and  4 at  the  points  7=36.5  and 
39.5  s. 

That  means  that  for  the  given  value  of  the  probe  nose 
radius  it  is  possible  to  duplicate  in  the  IPG-4  plasmatron 
all  three  conditions  (3.3),  (3.4)  for  the  one  trajectory 
point  only.  For  the  nose  radius  R^=0.8  m the  point  of  the 
complete  heat  transfer  simulation  is:  36.5  s, 

77=45.4  km,  Va=5524  m/s  (M„= 28.9);  the  corresponding 
subsonic  jet  parameters  are  the  following:  />/=5.8-10'2atm 
and  7i/=15.3  MJ/kg.  For  the  nose  radius  Tf^l.25  m the 
desired  point  is:  t2= 39.5  s,  77=41.8  km,  FK=5163  m/s 
(Mo=26.7);  the  corresponding  subsonic  jet  parameters 
are  the  following:  p/=7.8-10"2  atm  and  /z/=13.3  MJ/kg. 
We  see  that  both  points  tt  and  t2  are  close  to  the  point  t 
where  the  heat  flux  has  a maximum  value  and  ts<t  <t2. 
At  these  two  trajectory  points  the  stagnation  point  heat 
fluxes  to  the  probe  would  be  less  than  the  maximum 
values  of  the  heat  fluxes  but  these  differences  would  be 
within  1%  (see  Fig.  10).  So,  in  fact,  a remarkable 
opportunity  for  complete  simulation  of  the  heat  transfer 
to  the  Mars  probe  quite  close  to  the  most  important  point 
of  the  entry  into  Martian  atmosphere  by  using  the  IPG-4 
plasmatron  and  the  euromodel  arises. 

We  should  emphasize  once  more  that  partial  duplication 
of  only  two  parameters  ps  and  hs  along  a trajectory  in  the 
real  time  scale  would  be  not  acceptable.  In  this  way  in 
test  we  shall  have  excessive  heat  fluxes  if  t<f  and 
underestimated  heat  fluxes  if  t>f,  where  f=tt  if  Rg=0.8 
m and  f=t2  if  R^r  1-25  m.  Therefore,  the  most  correct 
test  technique  for  the  study  of  the  thermochemical 
resistance  of  the  TPM  for  the  vehicle  during  his  entry 
into  Martian  atmosphere  would  be  tests  at  constant 
values pf  and h"  during  15  s. 


8.  Prediction  of  Subsonic  Carbon  Dioxide  Test  for 
Mars  Pathfinder 

Let  us  consider  another  application  of  the  LHTS  concept 
to  a prediction  of  the  high-enthalpy  test  conditions  in  the 
IPG-4  plasmatron  for  precise  simulation  of  the  stagnation 
point  heating  for  the  Mars  Pathfinder  aeroshell  which  is  a 
70-deg  sphere  cone  with  a nose  radius  of  0.6625m  [ref. 
3].  In  accordance  with  ref.  3 the  maximum  of  the  heat 
flux  is  achieved  at  the  altitude  40.7  km  and  the  following 
free  stream  conditions  in  the  Martian  atmosphere 
(Cco2= 0-97,  and  CN2= 0.03):  FK=6590  m/s  and 

^=3 .23 TO'4  kg/m3.  The  IPG-4  operating  envelope  for 
the  subsonic  regime  with  carbon  dioxide  as  a working 
gas  contains  the  peak-heating  parts  of  the  Mars 
Pathfinder  [ref.  3]  trajectory  (Fig.  9). 

The  desired  test  conditions  are  determined  for  a subsonic 
dissociated  carbon  dioxide  flow  around  the  cylindrical 
model  of  50  mm  in  diameter  with  a flat  face  and  the 
rounded  edge  of  1 1 mm  in  radius  which  is  supposed  to 
be  exposed  for  testing  in  the  100-kW  IPG-4  plasmatron 
with  the  quartz  discharge  channel  of  80  mm  in  diameter 
[refs.  1 1-14], 

Now  we  can  recalculate  the  entry  parameters  specified 
above  to  plasmatron  test  conditions  by  using  formulae 
(3.3),  (3.4)  and  (7.2).  This  simple  technique  (when  the 

effective  radius  R*m  is  known)  gives  the  following  test 

conditions:  the  enthalpy  77e=21.73  MJ/kg,  the  static 
pressure  ps= 0.14  atm,  the  flow  velocity  V,= 1 84  m/s.  For 
the  considered  test  configuration  the  calculated 
dimensionless  parameters  in  equations  (4.1)  are:  zl=0.4, 
wle=0.39,  ae=2.10. 


For  these  flow  parameters  and  test  geometry  described 
above  the  stagnation  point  heat  transfer  rates  have  been 
calculated  as  a function  of  the  surface  temperature  Tw 
and  the  effective  catalytic  efficiency  yw  from  the 
numerical  solution  of  the  ID  boundary  layer  problem 
(4.1),  (4.2)  for  a 5-species  dissociated  carbon  dioxide 
mixture  ( C02  , 02 , CO,  O,  and  Q. 

The  next  assumptions  have  been  made:  1)  the  following 
reactions  are  running  in  the  mixture:  C02+M<->C0+ 

+0+M,  o2+m<^o+o+m,  co+m-^c+o+m,  co+o<-> 

oC+02„  C02+0<b>C0+02;,  2)  the  surface  catalytic 
recombination  of  the  CO  molecules  in  the  reaction 
C0+0^C02  and  O atoms  in  the  reaction  0+0->02  are 
the  reactions  of  the  first  order  with  equal  efficiencies  yw, 
and  the  C atoms  are  not  involved  in  surface  reactions. 
The  same  fourth-order-accurate  finite-difference  scheme 
was  exploited  for  the  numerical  solution  of  the  boundary 
layer  problem  (4.1),  (4.2). 

The  computed  stagnation  point  heat  flux  envelope  for  the 
predicted  subsonic  test  is  shown  in  Fig.  12.  The  upper 
border  of  this  envelope  (the  curve  7)  corresponds  to  the 
fully  catalytic  surface  (yw=  1),  the  lower  border 
(the  curve  7)  - to  the  noncatalytic  surface  (yw=  0).  The 
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solid  curves  2-6  correspond  to  the  constant  values  of  yw, 
=10"',  3 - 1 0"2,  10'2,  310'3,  10"3,  the  line  8 corresponds  to  the 
theoretical  minimum  of  the  heat  flux  from  the  frozen 
boundary  layer  to  the  noncatalytic  wall. 

The  heat  flux  envelope  is  limited  from  the  right  side  by 
the  curve  qw=zthoTw4,  where  z,h  is  the  total  hemispherical 
emissivity,  cr  is  the  Stefan-Boltzmann  constant.  The 
curves  9 and  10  correspond  to  radiative-equilibrium 
walls  with  z,h=  1 and  0.78. 


Fig.  12.  Heat  flux  envelope  for  the  IPG-4  plasmatron  subsonic 
regime  for  the  duplication  of  the  stagnation  point  heat  transfer 
to  the  Mars  Pathfinder  aeroshell  (Rn=0.6625  m)  at  the  peak- 
heating point  (h=40.7  km,  Vx-6.59  km/s)  by  using  an 
euromodel  (Rm=2.5-I0'2m). 

Now  it  is  very  easy  to  determine  the  heat  flux  and  the 
radiative-equilibrium  surface  temperature  for  the  given 
values  of  z,h  and  yw'-  qw  and  Tw  are  just  the  coordinates 
of  the  intersection  point  for  two  curves  - qw^qwiJw, 
const)  and  qiV=zthcjT w' . Thus,  in  the  predicted  test  at 
s,/f=0.78  for  the  fully  catalytic  wall  we  have 
qw=  142  W/cm2  and  7)^2380  K,  for  the  noncatalytic  wall 
- qvr 47  W/cm2  and  7V=  1800  K. 

For  the  heat  transfer  rates  at  the  stagnation  point  of  the 
Mars  Pathfinder  aeroshell  and  the  radiative-equilibrium 
wall  we  have  the  maximum  values  qw=\21  W/cm2  and 
7V=2315  K at  the  wall  condition  of  fully  recombined 
C02  and  qw= 42  W/cm2  and  7V=1755  K in  the  non- 
catalytic wall  case  [ref.  3]. 

So,  we  observe  quite  sufficient  agreement  between  the 
numerically  predicted  whole  heat  flux  ranges  for  the 
Mars  Pathfinder  aeroshell  at  the  trajectory  peak-heating 
point  and  for  test  with  an  euromodel  in  the  subsonic 
high-enthalpy  carbon  dioxide  jet,  which  could  be 


performed  in  the  IPG-4  plasmatron  with  free  stream 
parameters  determined  on  the  basis  of  the  LHTS  concept. 

It  is  important  that  the  enthalpy  and  pressure  values 
determined  above  belong  to  the  operating  envelope  of 
the  IPG-4  plasmatron  in  a subsonic  mode  (Fig.  9 ) and 
the  maximum  of  the  thermochemical  load  on  the 
stagnation  point  of  the  Pathfinder  aeroshell  could  be 
duplicated  precisely  enough. 

Recently  it  was  shown  in  ref.  24  that  quite  precise 
duplication  of  the  stagnation  point  heat  flux  to  the  Mars 
Probe  [ref.  17]  also  could  be  achieved  by  using  the  IPG- 
4 plasmatron  in  subsonic  regime  with  carbon  dioxide  as 
working  gas. 

This  analysis  was  performed  without  taking  into  account 
the  surface  ablation,  but  in  general  the  LHTS  concept 
validated  here  can  be  applicable  to  a heat  transfer 
problem  with  ablation  effects  if  the  influence  of  injection 
into  the  incoming  flow  is  not  very  strong. 

9.  The  Range  of  Applicability  of  the  LHTS  Concept 

The  derivation  of  conditions  (2.7)  - (2.9)  of  the  local  heat 
transfer  simulation  was  based  on  the  boundary  layer 
theory  and  the  assumption  of  local  thermodynamic 
equilibrium  in  the  subsonic  flow  (outside  the  boundary 
layer),  which  is  fairly  accurately  satisfied  for  molecular 
gases  (air,  nitrogen,  oxygen,  carbon  dioxide)  in  inductive 
plasmatrons  at  pressures  p>0.1  atm. 

It  was  found  in  ref.  25  that  the  boundary  layer  theory  is 
applicable  for  calculating  stagnation  point  heat  transfer 
in  subsonic  high  enthalpy  flows  at  Reynolds  numbers 

Res  = P gVsR*m  / > 30 , if  the  velocity  gradient  at  the 

stagnation  point  on  the  model  is  determined  with  making 
allowance  for  the  finite  thickness  of  the  boundary  layer. 
Practically,  for  the  models  used  in  the  heat  transfer  tests 
in  the  IPG  plasmatrons  30<Re5  <3 -102. 

Thus,  the  general  conditions  of  the  LHTS  concept 
formulated  above  are  applicable  for  subsonic  high 
enthalpy  molecular  gas  flows  in  plasmatrons  when 
/?>0.1  atm  and  Res  >30.  For  the  surfaces  with  a high 
catalycity  the  LHTS  concept  is  valid  in  the  full  range  of 
the  pressure,  for  the  surfaces  with  moderate  catalycity  the 
lower  border  of  the  pressure  range  should  be  estimated. 
For  the  accurate  computations  of  the  velocity  gradient  or 
effective  radius  for  subsonic  test  conditions  the 
numerical  solution  of  the  Navier-Stokes  equations  must 
be  used. 

The  accuracy  of  the  prediction  of  ground  test  parameters 
or  extrapolation  to  flight  could  be  improved  if  the 
velocity  gradient  for  flight  conditions  is  calculated  more 
accurately  directly  from  the  numerical  solution  of  the 
hypersonic  shock  layer  problem. 
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10.  Comments  to  Thermochemical  Simulation  in 
Supersonic  Tests 

Traditionally,  in  high  enthalpy  tests  practice,  especially 
in  tests  performed  by  using  arc-jet  facilities,  supersonic 
regimes  are  used  more  often,  then  subsonic  ones.  The 
LHTS  concept  gives  us  the  strict  inequalities  for  the 
similarity  parameters  % and  £ (see  (2.10)),  when 
supersonic  high  enthalpy  flow  should  be  performed  for 
the  providing  of  the  correct  duplication  of  the  stagnation 
point  heat  transfer: 


y.s~l 

'^  + 1 


<c<  1 (10.1) 


These  inequalities  mean  that  some  supersonic  test  is 
necessary,  if  the  nose  radius  and  the  model  radius  are 
comparable,  although  Rm<RN  : practically,  for  the 
simulation  of  the  hypersonic  heating  of  a vehicle  with  a 
small  nose  radius. 

For  the  case  of  duplication  of  the  heat  transfer  to  a nose 
cap  with  small  radius  the  simulation  conditions  for 
velocity  and  enthalpy  (2.7)  remain  valid,  but  the 
conditions  for  pressure  (2.8)  must  be  modified.  For 
supersonic  simulating  flow,  in  which  the  gas  between  the 
shock  wave  and  the  outer  edge  of  the  boundary  layer 
near  the  flow  axis  is  in  the  equilibrium  state,  this 
modification  consists  in  the  following.  In  the 
approximate  Poisson  adiabatic  equation  (2.5),  applied  to 
the  equilibrium  gas  moving  behind  the  shock,  the 
pressure  and  Mach  number  are  eliminated  by  means  of 
the  Rankine-Hugoniot  conditions  on  the  shock  wave,  and 
equations  (2.7).  As  a result  we  obtain 


PS  _ 

PW  \ + (2itS(\-zs)-\)?' 


(10.2) 


Where  y,  is  the  effective  specific  heat  ratio,  subscripts  S 
and  Sh  relate  to  the  parameters  in  the  plasmatron  free 
stream  and  behind  the  shock. 

The  validation  of  the  conditions  (2.7)  and  (10.2)  as  the 
basis  of  the  LHTS  concept  for  supersonic  tests  is  not  in 
fact  yet,  because  the  corresponding  procedure  requires 
the  computations  of  the  high  enthalpy  reacting 
supersonic  flows  in  thermal  and  chemical  nonequilibrium 
in  the  framework  of  the  full  Navier-Stokes  equations. 


11.  Discussion 

The  accuracy  of  the  prediction  for  the  thermochemical 
action  of  the  reacting  gas  on  a vehicle  surface  for  a 
hypersonic  flight  conditions,  based  on  the  LHTS 
concept,  depends  on  the  surface  catalycity  and  the 
displacement  of  gas  temperatures  from  the  equilibrium 


values  at  the  outer  edges  of  the  boundary  layers  on  a 
body  and  a model.  For  surfaces  with  high  and  moderate 
catalycity  the  satisfaction  of  the  conditions  (3.3),  (3.4)  or 
(3.6)  ensures  the  accurate  duplication  of  the  convective 
heat  fluxes  and  the  diffusive  fluxes  of  atoms  in  subsonic 
high-enthalpy  jets,  when  R„,«RN. 

For  the  real  TPM  with  catalytic  efficiency  yf(/»3  ■ 1 0'3  the 
above  mentioned  conditions  of  a local  heat  transfer 
simulation  can  guarantee  in  test  not  only  full-scale  heat 
transfer  rates,  but  real  nonequilibrium  chemistry  within 
boundary  layer  as  well.  If  the  states  of  the  dissociated  gas 
flow  at  the  edge  of  the  boundary  layer  in  the  test  or 
hypersonic  flow  conditions  are  significantly  non- 
equilibrium, ground  test  parameters,  predicted  by  using 
of  the  LHTS  concept  in  order  to  simulate  the  heat 
transfer  to  a noncatalytic  surface,  may  be  applied  as 
support  values. 

The  LHTS  concept  makes  an  analysis  of  heat  transfer 
parameters  along  trajectory  rather  simple  and  gives  clear 
algorithm  of  the  determination  of  the  trajectory  point  for 
which  stagnation  point  heat  transfer  could  be  accurately 
duplicated  without  taking  into  account  the  actual 
information  about  TPM  catalycity. 

As  we  have  seen,  the  computed  heat  flux  range  for 
simulation  of  the  stagnation  point  heat  transfer  rate  to  the 
Mars  Pathfinder  aeroshell  at  the  trajectory  peak-heating 
conditions  in  the  subsonic  high-enthalpy  carbon  dioxide 
flow  is  found  in  sufficient  agreement  with  the  study 
carried  out  through  the  full  viscous  shock  layer 
computations  in  ref.  3. 

Conclusion  & Outlook 

It  is  common  knowledge  now  that  modem  ground 
facilities  are  unable  to  produce  all  flight  conditions 
above  Mach  8.  At  the  same  time  the  quantitative  heat 
transfer  could  be  duplicated  in  plasmatron  quite  precisely 
at  least  for  a vehicle  stagnation  point.  The  LHTS  concept 
reveals  the  new  capabilities  in  the  planning  of  a high- 
enthalpy  experiment  and  the  new  approach  to  the 
extrapolation  from  ground  to  flight.  For  the  complete 
stagnation  point  thermo-chemical  simulation  the  triad  of 
parameters  - total  enthalpy-stagnation  pressure-velocity 
gradient  - must  be  duplicated  in  a high-enthalpy  test.  In 
this  kind  of  the  heat  transfer  simulation  the  pairs  of  the 
parameters  M-Re  and  pL-V  [ref.  26],  which  are  widely 
used  in  aerodynamics,  are  not  the  similarity  parameters 
in  the  stagnation  point  heat  transfer. 

The  conditions  of  the  hypersonic  flow  past  a blunt  body 
and  the  desired  conditions  of  the  free  stream  in  ground 
test  could  be  easily  linked  on  the  basis  of  the  LHTS 
concept  if  one  knows  or  can  calculate  an  effective  model 
radius  by  using  CFD  for  the  test  configuration.  For  the 
typical  entry  trajectory,  if  the  geometry  of  a model  is 
specified,  those  conditions  determine  only  one  trajectory 
point,  for  which  the  complete  local  simulation  could  be 
achieved.  If  the  trajectory  point  and  the  nose  radius  are 
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specified,  conditions  of  the  heat  transfer  duplication 
determine  the  test  conditions  and  the  effective  radius  of  a 
model. 

In  the  more  dramatic  noncatalytic  surface  case,  we  have 
to  solve  numerically  not  only  the  nonequilibrium  shock 
layer  problem  for  flight  conditions,  but  also  to  compute 
the  nonequilibrium  plasma  flow  within  a plasmatron 
discharge  channel  and  a subsonic  (or  supersonic) 
reacting  gas  flow  past  a model  for  the  prediction  of  the 
well-documented  test  conditions.  So,  CFD  modeling  is 
an  indispensable  tool  for  the  construction  of  the  bridge 
from  ground  test  to  flight  and  for  the  verification  of  the 
quality  of  flight  parameters  duplication  in  plasmatron 
tests.  The  interaction  between  ground  testing  and  CFD 
modeling  is  a genesis  for  real  gas  effects  duplication  and 
extrapolation  to  flight. 

Above  examples  of  the  LHTS  concept  applications 
clearly  shown  that  the  maximum  thermochemical  load  on 
a vehicle  surface  at  the  stagnation  point  and  TPM 
behavior  can  be  directly  duplicated  for  the  wide  range  of 
the  reentry  and  entry  conditions  in  the  Earth  and  the 
Martian  atmospheres  by  using  the  inductive  plasmatrons 
in  the  subsonic  regimes. 
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